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The structure of short- and medium-range order of silicate glass containing additions of phosphorus oxide is 
simulated on a computer. The phosphorus atoms can be in oxygen tetrahedra with different symmetry as well 
as in pentahedra and trigonal bipyramids. Calculations of the P - O bond lengths and the azimuthal angles in 
the phosphate glass-forming polyhedral are performed and a conclusion is drawn concerning the most likely 
structural position of phosphorus. The shortcomings of purely ionic models are noted. 
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One trend in the modern glass industry is a transition to 
high-calcium compositions with low magnesium oxide con¬ 
tent, making it possible to form glass articles rapidly. While 
having indisputable economic and technological advantages, 
such a transition creates problems due to an increase in the 
crystallization capacity of the melt and decrease of the chem¬ 
ical stability and mechanical strength of the articles. This 
problem can be solved by introducing small amounts of addi¬ 
tives, specifically, B 2 0 3 and P 2 0 5 [1,2]. In a previous work 
[3] we introduced P 2 O s into commercial silicate glass, i.e., 
doping the silicate matrix with phosphorus, which gives defi¬ 
nite technological advantages. 

The effect of phosphorus oxide additives on the short- 
and medium-range order of the silicate melts and glasses, 
i.e., at the nanoscale level, can be substantial and has not 
been adequately studied thus far. This is due to, first and 
foremost, the complicated character of the interaction of 
P 2 O s with the main glass former Si0 2 , as well as with A1 2 0 3 
and the network modifiers. Here phosphorus can become in¬ 
corporated into the silicate framework or form phase micro¬ 
nonuniformities, which is reflected differently on the proper¬ 
ties of glass. 

We studied the structural characteristics of the short- and 
medium-range order in multicomponent silicate glasses con¬ 
taining additions of phosphorus oxide. This investigation is 
based on computer simulation. 

The most commonly used models of glass are based on 
molecular dynamics and Monte Carlo methods. However, 
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the conventional schemes of these methods are too inaccu¬ 
rate, even for Si0 2 and B 2 0 3 based glasses, because the spe¬ 
cific natures of the bridge and terminal oxygen atoms are ne¬ 
glected. The possibility of modeling P 2 0 5 -based glass is 
even more problematic. This is because there are three bridge 
oxygen atoms in the [P0 4 ] tetrahedra and the fourth is bound 
by a double bond to the central phosphorus atom. Such a 
structure cannot be described by the universal P - O interac¬ 
tion potential. 

According to the information presented in [4], when 
modifying additives are introduced into phosphate glass the 
number of double bonds decreases. This can occur as a result 
of symmetrization of the [P0 4 ] tetrahedra. The MPDP 
method (modified method of neglecting differential overlap) 
used by B. S. Vorontsov [5] made it possible to obtain the 
computed values of the bond lengths P - O and P = O as well 
as the azimuthal angles, which are all in good agreement 
with the experimental data [6] and nonempirical calculations 

[7] , The change of the optimal geometry and the equalization 
of the phosphorus - oxygen bond are confirmed on the basis 
of these models. 

Modeling phosphate systems by means of molecular dy¬ 
namics on the basis of data from low-angle x-ray diffraction 
using the Born - Mayer potential and ion charge + 5 and - 2 

[8] does not make it possible to obtain a predominately tetra¬ 
hedral coordination of phosphorus. A substantial number of 
the P 5+ ions are at the center of five-vertex bipyramids, 
which are interconnected via common vertices, the P - O - P 
angles being close to the Si - O - Si angles in silica. 
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TABLE 1. 


Type of polyhedron 

Bond length, A 

Bond angle, ...° 

Tetrahedron with 

O 

II 

1 

O 

II 

1 

O 

1 

phosphoryl bond 

1.60282 

106.7640 


1 

o 

1 

hd 

1 

110.8760 


1.68206 

1 

O 

1 

1 

O 

1 


1.68264 

110.2510 


1.68677 

108.9580 

Tetrahedron with 

1.68764 

107.6550 

equivalent bonds 

1.68695 

109.6330 


1.68420 

111.8240 


1.69287 

108.1070 

Trigonal bipyramid 

1.98432 

82.3927 


1.83254 

100.7540 


1.95308 

69.1280 


1.72995 

123.3730 


1.80998 

70.4437 

Pentahedron 

Pyramid base 


1.75649 

85.0164 


1.75500 

81.3462 


1.75828 

101.9050 


1.75134 

86.4281 


Vertex 


1.74821 

106.2360 


- 

94.1046 


- 

87.8930 


- 

109.1110 


We have constructed a model of glasses in the system 
Si0 2 - P 2 0 5 - A1 2 0 3 - Na 2 0 using the Giperkhem-7 pro¬ 
gram in the ionic-bond variant, which makes it possible to 
obtain better matching with the previously established struc¬ 
tural characteristics. It should be noted that the applicability 
of the ionic model to phosphate glass cannot viewed as com¬ 
pletely justified. The Pauling electronegativities of the O and 
P atoms are 3.5 and 2.1, and for the difference of 1.4 the 
iconicity of the P - O bond is only 39%. This is appreciably 
lower than in silica (51%). 

Our model allows phosphorus to be in an oxygen envi¬ 
ronment of different symmetry. For phosphorus several dif¬ 
ferent variants of the short-range oxygen environment were 
considered: for coordination number 4 — a tetrahedron with 
a double (phosphoryl) bond and tetrahedron with equal 
bonds tied with a [A10 4 ] tetrahedron and for CN = 5 — 
pentahedron and trigonal biopyramid. 

Minimization of the geometry shifted the atoms from 
their initial positions. All optimized phosphate glass-forming 
polyhedral differ by the degree of departure from symmetry. 
The main parameters of the topological models of phosphate 
glass-forming polyhedral are presented in Table 1. 

It follows from the data in Table 1 that here the diffe¬ 
rence between the P = 0 (1.6 A) and P-0 (1.68 A) bond 


O 



Fig. 1. Structural fragments of silicate glass, containing quadruply 
coordinated phosphorus with double (a) and equivalent ( b ) bonds. 

lengths is appreciable, though this difference is much smaller 
than in pure P 2 0 5 (according to [4], 1.57 and 1.67 A, respec¬ 
tively). This attests to structural symmetrization of ortho¬ 
phosphate tetrahedra when modifiers are present in the glass. 
Figure 1 shows the structure of short-range order around an 
orthophosphate tetrahedron with double phosphoryl and with 
equivalent bonds. 

As result of the presence of uncompensated positive 
charge on the tetrahedron (Fig. lb ) its position is stabilized if 
A1 2 0 3 is present and ties can be formed with [A10 4 ] tetra¬ 
hedron, carrying a negative charge. Here the parameters of 
the bonds and the azimuthal angles different negligibly as 
compared with the parameters of the single bonds in a tetra¬ 
hedron with a double bond. 

It is evident that in the short-range order structure around 
pentacoordinated phosphorus, present in the form of a penta¬ 
hedron, the P - O bond lengths are shorter and this structural 
polyhedron is more symmetric than a trigonal bipyramid. For 
pentacoordinated phosphorus located at the center of a 
trigonal bipyramid it should be noted that the P - O bonds 
are systematically longer and the symmetry of the structural 
polyhedron is lower. 

Assessing as a whole the results of the construction of 
topological models describing the structural position of 
phosphorus in silicate glasses with high content of modifiers 
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Fig. 2. Structure of short-range order around pentacoordinated 
phosphorus in the form of a pentahedron ( a ) and trigonal bipyra¬ 
mid (b ). 

it is evident that such models can be realized, tough with dif¬ 
ferent degrees of likelihood. Since glass is amorphous and 
has a strongly nonequilibrium structure, which forms during 


cooling of viscoelastic melt, the structure of the nearest envi¬ 
ronment of phosphorus can differ substantially from similar 
structures in crystalline compounds. On the basis of energy 
considerations, structures with shorter bond lengths, i.e., te- 
trahedra, are more likely to be formed. It is less likely to find 
phosphorus in oxygen polyhedra, where its coordination 
number CN = 5. The formation of pentahedra is preferable. 

It should be noted that models employing ion charges + 5 
and - 2 cannot be accurate for such complex objects as 
multicomponent silicate glasses doped with phosphorus. It 
would be more promising to construct models in an ionic - 
covalent variant. The introduction of a mixed bond makes it 
possible to improve the correspondence between the empiri¬ 
cal data on the properties and the structural characteristics. 
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